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Abstract. Conventional X-ray methods use incoming plane waves and result in discrete diffraction 
patterns when scattered at crystals. Here we find, by a systematic method, incoming waveforms 
which exhibit discrete diffraction patterns when scattered at helical structures. As examples we 
present simulated diffraction patterns of carbon nanotubes and tobacco mosaic virus. 

The new incoming waveforms, which we call twisted waves due to their geometric shape, are 
found theoretically as closed-form solutions to Maxwell’s equations. The theory of the ensuing 
diffraction patterns is developed in detail. A twisted analogue of the Von Laue condition is seen 
to hold, with the peak locations encoding the symmetry and the helix parameters, and the peak 
intensities indicating the electronic structure in the nnit cell. 

If suitable twisted X-ray sources can in the future be realized experimentally, it appears from 
our mathematical results that they will provide a powerful tool for directly determining the detailed 
atomic structure of numerous biomolecules and nanostructures with helical symmetries. This would 
eliminate the need to crystallize those structures or their subunits. 
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1 Introduction 


This paper explores - at the level of modelling and simulation - the possiblity of novel X-ray 
methods for the determination of the detailed atomic structure of highly regular but not periodic 
molecules. The details are worked out for helical structures. These include carbon nanotubes, 
the necks and tails of viruses, and many of the common proteins (actin, collagen). The quest for 
novel methods is motivated by the fact that current X-ray methods, while hugely successful, have 
important shortcomings. A native helical assembly of proteins either has to be broken at the outset 
and the proteins crystallized, which is difficult and may lead to non-native forms; or one uses X-ray 
hber diffraction, which resolves only the axial but not the angular symmetry into sharp peaks. 

Roughly, our idea is the following. Conventional X-ray methods use incoming plane waves 

E(x,t) = B(x,t) = ;^(k x (1) 


and result (in the relevant regime of X-ray wavelength << sample diameter << distance of detector 
from sample, Fresnel number << 1) in the outgoing held 


Eoiii (x, t) 


const 

X - Xc 




[ dy, 

Jn 


( 2 ) 


with polarization vector n' = (I — -j^ ® |^)n and outgoing wavevector k'(x) = |k||^5^- Here p is 
the electron density of the illuminated sample and Xc denotes a typical point in the sample. 

The emergence of the Fourier transform in (2), and its amazing properties regarding construc¬ 
tive/destructive interference, underlie the power of X-ray methods for periodic structures. One can 
see from (2) and its derivation that the Fourier integral kernel -O y jg directly arising from 
the assumption of a plane-wave source (1). Other sources would give other kernels. This suggests 
the following line of research: design the ineoming radiation (as a solution of Maxwell’s eguations) 
such that the kernel interacts with highly symmetric hut non-crystalline structures with the same 
dramatic properties of constructive/destructive interference as occurs in the periodic case. 
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This design problem for the incoming waves can be formalized into the following mathematical 
problem: find time-harmonic solutions to Maxwell’s equations which are simultaneous eigenfunc¬ 
tions of a continuous extension of the generating symmetry group of the structure. Why this is a 
good formalization is a long story, told in Section 

For discrete translation groups, which are the generating symmetries of crystals, we show that the 
design problem is solved precisely by the plane waves used in classical X-ray methods. This is a 
new characterization of plane waves. It explains why plane waves are right for crystals. 

For helical symmetry groups, the design problem can also be completely solved. The ensuing 
family of incoming waves is 




cos ip — sin ip 0 
sin p cos p 0 

0 0 1 



Ja+iinr) 

Ja-lM 
JaM 


(3) 


where (a,/?, 7 ) G Z x M x ( 0 , 00 ) is a parameter vector analogous to the wavevector k in Q, 
n G is a polarization vector which must satisfy ( 0 , 7 ,/ 3 ) ■ n = 0 , and the frequency u is given 
by ca = c|( 0 , 7 ,/5)|. The cartesian vector ( 0 , 7 ,/ 9 ) has a simple physical meaning which will emerge 
The Ja are Bessel functions, (r, p, z) are cylindrical coordinates with respect to the 


in Section 14 


helical axis, and Eg is the cartesian field vector, with the third component corresponding to the 
axial direction. We call the electric fields (|^ twisted waves. Figure shows the twisted wave with 
parameter vector (a, /?, 7 ) = (5, 3,1) and polarization vector n = (1, 0, 0). 



Figure 1: A twisted wave with angular, axial and radial wavenumber (q:,/ 3 , 7 ) = (5,3,1). The plot 
shows the real part of the first component of the wave. Red: positive values; Blue: negative values. 
Note the helical shape of the level sets restricted to a co-axial cylinder, and the Bessel pattern 
perpendicular to the axis. 

Thus twisted waves consist of four factors: a scalar plane wave on the cylinder; a rotation matrix 
which rotates the field direction along with the base point; a somewhat mysterious polarization ten- 
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sor which depends on the polarization vector; and a vector of three Bessel functions of neighbouring 
order. 

We remark that, while the axial component is just a scalar cylindrical harmonic, the twisted 
wave ([^ is not a cylindrical vector harmonic, except for special choices of nj^ 


The helical shape of the level sets in Figure suggests that for suitable values of the angular and 
axial wavenumber, twisted waves can induce resonant electronic oscillations of every single molecule 
in a structure with helical architecture. Thus one can hope for diffraction intensities which strongly 
depend on the twisted wave parameters. 

At least in axial direction, the radiation scattered by a helical structure indeed exhibits sharp 
discrete peaks with respect to the radiation parameters a and /3. More precisely: the signal of 
a helical structure in axial direction vanishes unless the angular/axial wavenumbers of the twisted 
wave minus the axial wavenumber of the outgoing wave belong to the reciprocal helical lattice shifted 
left or right by precisely one angular wavenumber. This is an analogue of the Von Lane condition, 
but waves and structure are curved. The shifts come from the fact that the polarization direction 
of a twisted wave rotates along with the base point. Details are given in Section 11 Moreover, as in 


X-ray crystallography, the unit cell electron density can be recovered, up to a scalar phase problem, 
from the peak intensities. A further attractive feature is that the outgoing signal is invariant under 
axial translations and rotations of the structure^] 

These results suggest a hypothetical set-up of structure analysis with twisted X-rays. Send a 
twisted wave towards a co-axial helical structure. Use a detector further along the axis to record 
the diffracted intensities as a function of the incoming radiation parameters. Solve a scalar phase 
problem to infer the electron density. See Figure]^ 
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Figure 2: Hypothetical set-up of structure analysis with twisted X-rays. 

Twisted X-ray waves, (|^, are a theoretical proposal. Intriguingly, similar waveforms have been 
experimentally realized, such as optical higher-order Bessel beams |AD00j and photons and beams 
carrying an angular wave factor 6*““^ |AB92| IHD921 IMTTOTj . Clarifying the precise relationship 
to twisted waves would require an understanding of the full electromagnetic held of the produced 
modes. A related recent development is the creation of electron vortex beams |UT101 IVTSIO] . 

^The latter fields, introduced by Hansen |Ha34j . are defined as curUa-^) and curl curl (ab), where b is a scalar 
cylindrical harmonic and a S is a fixed “pilot vector”. This construction provides useful basis functions for the 
time-harmonic Maxwell equations, but it does not capture the geometric behaviour of the polarization of twisted 
waves. 

^By comparison, the signal poduced by fiber diffraction as described by the Cochran-Crick-Vand formula |CCV52] 
is not invariant under axial rotations, causing well-known difficulties in the interpretation of fiber diffraction images. 
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The plan of the paper is as follows. Sections 2-4 extend the standard scalar Fonrier transform 
model for plane-wave diffraction to a vector-valned electromagnetic model needed to treat general 
incoming waveforms. Sections 6-9 formulate the radiation design problem, and derive plane waves 
and twisted waves from it. Sections 10-12 develop the theory of diffraction patterns of twisted 
X-rays. Finally, in Section 13 we present simulated diffraction patterns of a carbon nanotube and 
of tobacco mosaic virus. 


2 The time-harmonic Maxwell equations 

Incoming electromagnetic waves will be sought as solutions to Maxwell’s equations in vacuum, 

1 9E 9B 

— — =curlB, divE = 0, —- =-curlE, divB = 0, (4) 

ot ot 

which are time-harmonic, that is to say 

E(x, t) = Eo(x)e-*‘^*, B(x, t) = Bo(x)e-'-* (5) 

for some a; > 0. Here E and B are the electric and magnetic helds and c is the speed of light. SI 
units are used throughout. The helds are dehned on x [0, cxd) and take values in C^. 

The ansatz ([^ reduces Maxwell’s equations Q to 

(^2 

AEq =-yEo, div Eq = 0, (6) 

Bo = --curlEo. (7) 

UJ 

We are interested in bounded solutions to 0 - We note that bounded weak solutions, that is to 
say vector helds in the function space C^) which solve 0 in the sense of distributions, are 

automatically bounded inhnitely diherentiable vector helds which satisfy 0 classically. 

A basic example of time-harmonic radiation, and the one used in classical X-ray crystallography, 
are plane waves 

Eo(x) = Bo(x) = -(ko x 11 ) 6 **^“'^ (8) 

OJ 

which solve ([^-([^ if the wavevector ko G and the polarization vector n G satisfy 

|ko| = -, ko • n = 0. (9) 

c 

Note that by the hrst of these relations, the time frequency a; in (|^ hxes the spatial wavelength 
A = 27r/|ko| of plane-wave radiation as 

A = 2ti-. (10) 

UJ 

The more general form ([^ is important when scattering from non-crystalline structures is under 
consideration. In fact the notion of wavelength and the relation remain meaningful in this case. 
This is because the Fourier transform Eo(k) = e“*^ ’’‘Eo(x)(ix is, due to eq. (|^, supported on 

the wavevector sphere |k| = cu/c, and hence any solution Eo can be thought of as a superposition 
of plane waves with wave vectors of equal magnitude uj/c. 

In X-ray crystallography, due to the goal of atomic resolution, only A’s comparable to interatomic 
distances, i.e. below a few Angstrom, are useful. Electromagnetic waves in this regime are known 
as hard X-rays. 
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3 Electromagnetic model for the diffracted radiation 

The standard scalar Fourier transform model (or oscillator model) of X-ray diffraction patterns 
which underlies the crystallography and biocrystallography literature is insufficient for our pur¬ 
poses. This is because it neglects the polarization of the incoming wave from the outset. However, 
fluctuating polarization directions turn out to be a necessary feature of any incoming field which 
can induce resonant electronic oscillations in a helical (or any other symmetric but non-crystalline) 
structure. As a consequence we need a full electromagnetic (vector-valued) extension of the standard 
model. 

Before presenting this extension, we briefly recall the standard scalar model, for comparison and 
background. 

Oscillator model. As described in many texts (see e.g. |AM76] ). one starts from a simplified 
picture which involves only scalar waves instead of electromagnetic fields and bypasses Maxwell’s 
equations. A scalar incoming plane wave hitting an electron at position y is assumed to 

induce a circular wave with same wavelength and same phase at the point y, 

const |x-y|-<^t)giko-y_ ( 11 ) 

Suppose that the incoming wave travels at the speed of light, |ko| = co/c; let Xc be a typical point 
in the sample; approximate the phase |ko||x — y| of the outgoing circular wave at a point x in the 
far field by k'(x) ■ (x — y), with direction-dependent outgoing wavevector 


k'(x) 


U X — Xc 

C |x — Xc\ 


( 12 ) 


and take a superposition of electrons with density p(y). This yields the scalar outgoing field 


^Pouti^,t) = const [ p(y)e-*(*^'W-*‘“)'ydy. 


Note that the integral appearing above, 

/(k'-k„)= [ e-‘l‘‘'-'‘”>>-p(y)<iy = p(k'-k„), 


(13) 


(14) 


is just the Fourier transform of the electron density, evaluated at the difference k' — k of outgoing 
and incoming wavevector. This integral encodes the density of the structure one wants to image, 
and is known as the structure factor or form factor or scattering factor. The squared amplitude 
(x, f)p of the scalar field (13) is then taken as a model for the measured intensity of the 

one obtains 

^o)P. (15) 


^(x) = _ 

outgoing electromagnetic radiation. From 

/(x;ko) = |V’o«t(x,f)P = |p(k'(x) 


Thus the “intensity of diffracted radiation” at the detector position x is the absolute value squared 
of the Fourier transform of p, evaluated at the difference k'(x) — ko between outgoing and incoming 
wavevector. To emphasize the dependence of I on both the observation point (which determines the 
outgoing wavevector) and the parameters of the incoming wave we have used the notation /(x; ko). 
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In particular, the data gathered from a sufficiently large set of observation points x and incoming 
wave parameters kg delivers the abstract data set 


/(k) = |Mk)p, ke 


(16) 


and the X-ray interpretation problem, i.e. the task of inferring atomic structure, p, from X-ray 
diffraction data, I, reduces to the phase problem for the Fourier transform. Note that the data 


set (16) is a function on the dual (wavevector) space of the physical space on which the density 


is dehned. The abstract “diffraction intensity” or “diffraction spectrum” or “diffraction measure” 


(16) constitutes the starting point of previous mathematical work on X-ray diffraction of crystals 


|St941 IFrO?] and quasicrystals |Ho95( IBMn4( IBGn8( IBG13j . 


The semi-empirical oscillator model (12), (13), (15) is sufficient for a great deal of X-ray physics. 


In particular, it accounts for the fluctuating phases of incoming plane-wave radiation 
Electromagnetic model. A proper vector-valued expression for the outgoing electric held which 


improves the model (13) and can deal with huctuating incoming polarization directions can be 


obtained as follows. Start from the well known Lienard-Wiechert helds (see e.g. |Gr99j i of a 
moving point charge driven by any given incoming held; pass to a non-relativistic (weak-held) limit; 
take a superposition of charges; and make a far held approximation corresponding to the regime X- 
ray wavelength << sample diameter << distance of observation point from sample, Fresnel number 
<< 1. The details can be found in our companion paper |F,1,115] and we only give the resulting 
expression for the dihracted electromagnetic held: 


(i - la ® la) F 

X Eont(x,t), 


^out 


(17) 


with direction-dependent outgoing wavevector familiar from the oscillator model. 


k'(x) = 


Ul X 


C X 


(18) 


Here Eo(|/)e“*“'* is the incoming electric held, a solution to the time-harmonic Maxwell equation 
(|^, p : —)■ M is the electron density of the sample, x^ is a typical point in the sample, and 

Cee is a universal constant depending, among other things, on the charge and mass of the electron. 
Moreover I is the 3x3 identity matrix and a® b denotes the 3x3 matrix A with entries Aij = aibj, 
where a and b are any two vectors in 

If the incoming electromagnetic held is replaced by its real part, as it properly should to model 
physical incoming X-rays, the dihracted radiation is given by the real part of 0- 

According to the model (17)-(p(^, the outgoing wavevectors k' have the same length as the 


incoming wavevectors k (recall from the previous section that any solution Eq to eq. (§ can be 
viewed as a superposition of plane waves with wavevectors k of equal length |k| = u/c), 


Ik'l = Ikl. 


(19) 


This relation has the important microscopic physical interpretation that the photon energy is con¬ 
served in the scattering. Note that the energy of a photon with wavevector k is E = |p|c, where 


p = hk is the photon momentum. Thus the model (17) corresponds to elastic or Thomson scatter¬ 
ing. 
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In the special case of incoming plane-wave radiation 
field reduces to 


-(§. eq. ( [T7| ) for the electromagnetic 




gi(k'(x)-x—(iip 


-n X 


X — Xr 


p(y) e 


-i(k'(x)-ko) 


■^dy, Bo„i(x,t) = -k'(x) x Eo„t(x,t), 

ijJ 


with outgoing direction-dependent polarization vector 

k'(x) ^ 


n X = 


I 


Ikhx) 


k^(x) 

Ikhx) 


n. 


( 20 ) 


( 21 ) 


In X-ray experiments the electromagnetic helds E and B cannot be measured directly; a detector 
only records the scalar held intensity. The latter is dehned (see e.g. |Gr99j ) as the time average of 
the absolute value of the Poynting vector, and has the physical dimension of power transferred per 
unit area, i.e. energy transferred per unit area per unit time. In formulae, 


/(x) = lim — [ |S(x, f)| df, S(x, f) = —E(x, f) x B(x, f). 
T^oo T Jo Ho 


( 22 ) 


Here and below, Ho and £o are the magnetic and electric constants, which are related to the speed 
of light by the equation EoHo = l/c^, and S denotes the Poynting vector. Note that the latter is 
the hux vector for the energy density e in Maxwell’s equation, that is to say real-valued solutions 
to (|^ satisfy 


d 

—e + divS = 0, where e(x, t) = - ( £o|E|^ H-|B 

dt 2 \ Ho 


(23) 


Here the restriction to physical, real-valued solutions is essential, and we do not recommend that 
the above dehnition of the Poynting vector be applied to complex helds. In Appendix 1 we evaluate 
expression (22) for the physical dihracted radiation, given by the real part of the helds Eo„t and 
Bout in eq. (17). The result, i.e. the intensity of the real part of the elecromagnetic held 0 at 


the observation point x, is 


-f(X = |Eooi(x, i)P. 


(24) 


where Bout is the complex electric held in (17) (note that its absolute value is independent of t). The 
at hrst sight mysterious appearence of the complex held amplitude comes from the time averaging 
in (22), as is clear from the derivation in Appendix 1. 


In the plane-wave case, the intensity can be read oh from eqs. (20) and (24). Noting that n' is the 
projection of n onto the orthogonal subspace of k', and assuming for the moment that n is real up 
to a phase factor (i.e., the wave is linearly polarized), we have 


and hence 


J(x;ko) = 


ceo 


X — X, 


n'l = sin (Z(n, k')) |n|, 

-—J siii^(Z(n,k'(x))) |np|p(k'(x) - ko)p. 


(26) 


(26) 


(Here we have used the same notation as in (15) to emphasize the dependence on the incoming 
wavevector kg.) If n is complex, the sin^ factor has to be replaced by |n'(x)p, with n' given by (21). 
In particular, the electromagnetic model (|T7|)-(@ specialized to plane waves recovers in a natural 




















way the well known fact - missed by the oscillator model - that scattering along the polarization 
direction of the incoming X-ray beam, corresponding to k'(x) being parallel to n, is completely 
suppressed. In the special case of a single electron, p = Sq, expression (26) agrees with what is 
typically assumed as a starting point in X-ray books (see e.g. |AMllj Section 1.2). Note also 
that the last factor in (26) is precisely the intensity of the semi-empirical oscillator model, i.e. the 
absolute value squared of the “structure factor” (14). 


Although the arguments in |FJJ15j leading to (20)-(21) are purely classical, the result agrees (up 
to the decay factor l/|x—x^, and granting the identihcation of macroscopic outgoing directions with 
microscopic photon wavevectors) with the result derived in the beautiful tutorial article |Sa09j from 
perturbative nonrelativistic quantum electrodynamics (NRQED). In particular, NRQED conhrms 
equations (18) and (21). 

In the sequel, we shall work with the electromagnetic model 


4 Generalized structure factor 


From the point of view of contemplating novel forms of incoming radiation, the most important 
outcome from the electromagnetic model ([T^-([T^ is a prescription for how the “structure factor” 
(14) needs to be modihed when the incoming X-rays are not plane waves. The corresponding term 


in (17) is the integral term 


fE„(k') = / Eo(y)p(y)e ’’‘’ ‘'dy 


(27) 


We call this term the generalized structure factor With respect to the incoming electric held Eq. The 
latter may be any bounded solution to the time-harmonic Maxwell equation ([^. Expression (27) 
is a vector-valued generalization of the familiar structure factor 


/(k' - ko) = / p(y) e 


-i{k'-ko)-: 


(eq. (14), Section 


(28) 


for classical (plane-wave) diffraction, and determines the amount of constructive or destructive 


interference of the outgoing waves in the direction k'. Note that the interference in (27) depends 


not just on the incoming and outgoing phases, but also on the incoming polarization directions. 

In the plane-wave case ([^, that is to say Eo(y) = the exponential factor from the incom¬ 


ing wave nicely combines with the exponential factor from (27), and the generalized structure factor 


reduces to just the classical structure factor (14), multiplied by the constant incoming polarization 
vector: 


fEo(k')=n / p(y)e 


-i{k'-ko) 


J-dy = n/(k'- k„ 


(29) 


It will be useful to interpret (27) mathematically as an integral transform that depends on the 


incoming held Eq and maps the electronic charge density p to a vector held on the dual (wavevector) 
space. It generalizes the Fourier transform which arises in the oscillator model, and might be called 
the radiation transform of p with respect to Eq. This transform is dehned by 


(77-EoP)(k) = / Eo(y)p(y)e •''■’'(iy (k e E®). 


(30) 


Here, as before, Eq may be any bounded solution to the time-harmonic Maxwell equation (|^. Eq. 
(30) naturally establishes 77 eo as a linear map from L^(]R^) to the space of bounded continuous 
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vector fields on the dual (wavevector) space Like the Fourier transform, it can be extended 
to tempered distributions, and maps these to vector-valued tempered distributions. This extension 
will be useful when discussing diffraction in the idealized but important case of inhnite helices, 
nanotubes and crystals. 


5 Classical Von Lane condition 


Given a complete mathematical model of the outgoing electromagnetic held such as , we can cast 
the Von Lane condition of classical X-ray crystallography, hrst discovered by Friedrich, Knipping, 
and Von Lane jFKL12j . in the form of a mathematical theorem. Let £ be a Bravais lattice in 
i.e. a set of the form 

C = AT? for some invertible 3x3 matrix A. (31) 

Let p be any smooth ^-periodic function on that is to say 

p(x -|- a) = p(x) for all x G a G £. 


By choosing a suitable partition of unity, any such p can be written in the form 

for some smooth, rapidly decaying function p belonging to the Schwartz space 5(M^). In the sequel, 
we will use the following notation for delta functions common in the mathematical literature: if xq 
is a point in and S a countable set of points in we write 


i5x«(x) i5s(x) =^i(x-a). 


(33) 


agiS 


In the Von Lane condition, the reciprocal lattice of C will naturally emerge. The latter is dehned 
by 

£' = {k G : k ■ a G 27rZ for all a G £}, (34) 


and is given in the case of (31) explicitly by 

C = 2'kA-^T}, 


(35) 


where A~'^ denotes the transpose of the inverse of the matrix A. The Von Lane condition can now 
be stated as follows. 

For simplicity, we only consider the square root of the intensity, i.e., up to trivial constants, the 
absolute value of the electric held (20). This is convenient because the held itself, unlike its square, 
is a well-dehned mathematical object for inhnitely extended systems such as (32) or p = 5£, by 
interpreting it as a distributional Fourier transform. The powerful machinery of Fourier analysis 
then allows to mathematically understand in a quick way the phenomenon of discrete dihraction 
patterns. 


Theorem 5.1 (Von Lane condition). Let C be the Bravais lattice (31), and /ef p : —)■ M be 
any C-periodic density, represented in the form (32). Assume that the incoming electric field is a 
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plane wave, Eo(y) = with polarization vector n G and wavevector ko G satisfying 

0 . and the outgoing radiation is given by the electromagnetic model with corresponding 

intensity (26). Assume moreover without loss of generality that = 0. Then 


(/(x;ko))^/^ = Co 


= Co 


det A\ 

(ML 

det A\ 


|x| 

1 


n — 


n — 


■ n)|^| 1^ (k'(x) - ko)| <5,:. (k'(x) - ko) 
|k'(x)|| l<^(a')l<^a'(k'(x) - ko), 


k'W . k'(x) 


.|k'(a;)| 


(36) 

(37) 


a'eC' 


where k'(x) = and Cq = In particular, the outgoing signal is zero unless the 

difference between outgoing and incoming wavevector, k'(x) — ko, is a reciprocal lattice vector. 


For the simple diffraction model (16), analogous results for were presented in |St941 [Ff-QT] . and 
a rigorous treatment of the intensity I renormalized by volume was given in |Ho95] . We remark 


that the right hand side of (36) makes rigorous sense (as a locally bounded measure) if considered 
as a function of the incoming wavevector ko at hxed observation point x, as implicitly assumed in 


the simple model (16), but not if considered as a function of x at hxed ko. 


Proof The central term in the outgoing held (20) is the structure factor fEo(k') = np(k' — ko). 


which is well-dehned as a distributional Fourier transform. It is convenient to re-write formula (32) 
as the convolution 

p = ip*6c, 

where (/ * (y')(x) = /(x — y)g{y) dy. By the generalized Poisson summation formula (see e.g. 

[Fr07j h the Fourier transform of dc is 


5/: = 


(2vr)^ ^ 

det ^1 


(38) 


The Fourier calculus rule f * g = f g now gives 

(27r^3 


P = 


det A\ 


ip6i 


/(k'(x) - ko) = 


(2x) 


det A\ 


{(pdc) (k'(x) - ko) 


Since the expression for the structure factor is, by inspection, a locally bounded measure, its absolute 
value is well dehned (see e.g. PISI), and corresponds in the above case to replacing the factor if 


by its absolute value. The result now follows immediately from (26). 


6 The euclidean group of isometries 

A key role in the following is played by the euclidean group of isometries 

E{3) = {(R|c) : R G 0(3), c G M^}, 

where 0(3) = {R G : R^R = 1} is the orthogonal group. Elements of E{3) act on points in 

three-dimensional euclidean space as 


(R|c)x = Rx -I- c. 
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The natural induced action on vector fields E : —)■ is 


((R|c)E) (x) = RE (R"^(x - c)) . 


(39) 


Note that R —c) is the inverse element (R|c) 
equations. 

Important continuous subgroups are 


The action (39) preserves solutions to Maxwell’s 


T={(I|c):cgM^}, xi-)-x + c (translation group) 


(40) 


and 


T-Le = {(Rke) : R G 50(3), Re = e, r G M}, x i—>• Rx + re (helical group with axis e), (41) 

where e is a given unit vector in Here 50(3) denotes the special orthogonal group, i.e. the 

elements of 0(3) with determinant 1. As we shall see shortly, the former subgroup is associated 
with classical plane-wave X-ray methods for crystals, and the latter with twisted X-rays. 


7 The design equations 

We now have a second look at plane waves. Why are they the right radiation to use for the analysis 
of crystals? 

Every researcher interested in X-ray diffraction is familiar with the Von Lane condition: the 
outgoing signal of a plane wave scattered at a crystal concentrates on a discrete set of outgoing 
directions, while in between it is suppressed by destructive interference. But suppose plane waves 
were not given to us a priori, as the radiation emitted by conventional X-ray tubes. Would we know 
how to come up with them if our goal was to achieve a nice diffraction pattern? What, exactly, is 
the “connection” between a particular family of solutions to Maxwell’s equations on the one hand, 
and point sets with crystalline order on the other? 

The connection is that plane waves and crystals have matching symmetries. By this we do 
not mean that they have the same symmetries. Plane waves have a larger, “continuous” family of 
symmetries. We first explain this informally, then make it precise in group-theoretical language, 
then generalize beyond crystals. 

Start with a crystal, i.e. a structure with atomic positions 

S = {xg^^ -|- a : a G £, i/ = 1,.., M}, C = {Wi -|- j \2 + : i, j, fc G Z} = A (42) 

where Xq^\ ••;Xg™'^ are the positions of finitely many reference atoms, vi, V 2 , V 3 are linearly inde¬ 
pendent vectors in and A is the matrix with columns given by these vectors. Thus the atom 
positions are obtained by translating finitely many reference atoms by each element of the lattice C. 
Mathematically, this means that the crystal is the “orbit” of a finite set of points under a discrete 
group C of translations, and in particular that each element of £ is a symmetry of the crystal, i.e. 
maps it to itself. The discrete translation group C is the generating symmetry group of the crystal. 
We now look at plane waves, 

Eo(y) = n G C^ ko G ko ■ n = 0. (43) 
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Plane waves also have a translation symmetry. Their valnes at different positions jnst differ by 
phase factors, 


Eo(xo + a) = (phase factor depending on a) Eo(xo) for all a G 


(44) 


with 


(phase factor in (44)) = e 


_ zko-a 


(46) 


Mathematically, eq. (44) means that the wave, a vector held on is an “eigenfnnction” of the 
operator which translates a vector held by a vector a G T_a : Eq e- >■ Eo(- + a). These operators 
form a continnous gronp which describes the action of the continnous translation gronp T on 
vector helds (as already discnssed in the previons section). Hence plane waves are simultaneous 
eigenfunctions of the continuous translation group T. 

We now come to the interaction between wave and strnctnre. The interaction occurs via the 


generalized structure factor (27) in the dihracted radiation held. Assume for simplicity that the 


density is a sum of delta functions at the atom positions, p{y) = X)ae£ ^xo+a(l/)- Then this factor is 


fEo(k') = / Eo(y)p(y)e = ^Eo(xo + a) 

ae£ 


—ik'a 


(46) 


The eigenfunction property (44) of the waves holds in particular for the crystalline translations 
a G £, and so it follows that 


fEo(k') = 


(phase factor depending on a) e 


—ik'a 


Eo(xo). 


(47) 


Thus the dihraction behaviour is reduced to a phase factor sum. The phase factors come from 
the symmetry of the wave. The points where they are evaluated come from the symmetry of the 
structure. And the sum “behaves nicely”: it interferes constructively when g-hk-ko) a _ ^ 
a G £, i.e. when k' — ko belongs to the reciprocal lattice C', and destructively otherwise. 

To get constructive interference, it would be enough if the wave had just the same symmetry as 
the crystal, i.e. if (44)-(45) was only true for a’s in the discrete translation group But to get 
destructive interference when the incoming radiation parameter ko or the lattice parameters A are 


tuned off resonance, one needs (44)-(45) for all a. 


In summary, the discrete diffraction patterns of classical X-ray crystallography can be traced to 
the fact that crystals and plane waves have matching symmetries. Constructive interference comes 
from the fact that plane waves share the symmetry of crystals. Destructive interference comes from 
the fact that plane waves have a larger, continuous symmetry group. 

Everything so far is just an abstract rationalization of a very well known phenomenon. But can it 
be generalized? 

The key is to realize that the foundation on which X-ray crystallography is built, the complex 


exponential form (43), can actually be derived from the innocent looking eigenvalue equation (44) 


^This condition has infinite-dimensionally many solutions: it just means that the electric field multiplied by a 
complex exponential, -ko) x^ shares the periodicity of the crystal. In the context of electron wavefunctions 

instead of electric fields, such waves are known as Bloch waves. 


13 













If we combine two translations, x i—>■ (x + a) i—>■ (x + a) + b, we can either apply (50) to the whole 
translation, or separately to the two translations by a and b, and so 


(phase factor at a + b) = (phase factor at a) ■ (phase factor at b). 


(48) 


Mathematically, this means that the phase factor is a group homomorphism from the translation 
group to the multiplicative group C\{0}. And the only such group homomorphisms are the scalar 
complex exponentials (45)! Eq. (45) together with (44) implies ([^ (with polarization vector 


n = Eo(xo)), up to the orthogonality constraint on n and ko- The latter follows from Maxwell’s 
equations, (|^. So plane waves are those solutions to (44) which satisfy Maxwell. (See Theorem 8.1 
for a precise statement.) 

Thus we have a path from crystals to plane waves. 


crystal —> generating symmetry group 

—> continuous extension of symmetry group 

—)■ eigenfunctions of continuous symmetry group which solve Maxwell. 


(49) 


This path can be generalized. We can start from any structure which is generated by a symmetry 
group, but the group does not need to consist of translations. Any discrete subgroup of the euclidean 
group E{2)) is hne. Many interesting structures in biology and nanotechnology have this form, 
with the generating subgroup often being a discrete helical group. The continuous extension of 
the group is then the helical group "He described in eq. (41). Solving the resulting combined 


eigenfunction/Maxwell problem will give a different family of incoming waves. 

The non-crystalline but highly symmetric structures which are generated by some discrete sub¬ 
group of E{3) form an interesting class. This class was recently introduced and studied by one 
of us |Ja06j . and has been named objective structures. Like crystals, objective structures can be 
completely classihed |DE,Tj . 


The remarkable constructive/destructive interference properties of the structure factor (46) sur¬ 
vive when structure and radiation are governed by non-translational symmetry groups. Here we can 
rely on the fact that mathematicians such as A. Weil |We64j realized a long time ago - albeit with 
a completely different motivation, number theory rather than molecular biology - that the Poisson 


summation formula (38) has far-reaching generalizations to sums of delta functions over discrete 


subgroups of continuous groups. For details we refer to Section ^ 


To conclude this section, we formulate the design problem suggested by the path (49) precisely. 


as a system of equations. As just discussed, crystals can be replaced by structures generated by a 
symmetry group which may contain rotations. 

Definition 7.1 (Design eqnations). Let G be a closed subgroup of the euclidean group E{3) 
- the “desired symmetry of the radiation”. (Typically, G is obtained as a continuous extension 
G D Go of a discrete subgroup Gq of E{3) - the “generating symmetry of a structure”.) A vector 
field Eq : —)• solves the design eguations for G if: 

(^)Eo is a joint eigenfunction of all group elements g G G; that is to say, 


gEo = x{g) Eq for some x{g) ^ C and all g E G. 


(50) 


Here the left hand side is given by eg. (39) which describes the action of the euclidean group E{3) 
on vector fields. 
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(ll) Eo is a bounded solution to the time-harmonic Maxwell equations 


ui 


AEo = —^Eo for some uj > 0, 
div En = 0. 


(61) 

(52) 


Note that the transformed held ^fEo (left hand side in (i)) reduces to the left hand side of (44) when 
the group element g is given by the translation (I| — a). 

In the next two sections, we will solve the design equations for some interesting examples. We 


will make use of the fact that as for (44), the eigenvalue xih) a function of g must be a group 
homomorphism: 

Lemma 7.1 (Character lemma). //Eq is any solution to the design equations for the group G 


which is not identically zero, then the function x ■ G ^ C in Def. |y. i| is a bounded continuous 
group homomorphism from G to the multiplicative group C\{0}. 


Proof The argument is the same as that leading to eq. (|4^. Eq. 


shows that for gi, g 2 ^ G, 


x(fi'ifi'2)Eo — (5'ifi'2)Eo — (7i((72Eo) — x{9i)92^o — x(5'i)x(fi'2)Eo. 


(53) 


Evaluation at a point x where Eq is not zero shows that x is a group homomorphism. Boundedness 
and continuity of x are a direct consequence of the same properties for Eq. 

In the case when G is abelian, the bounded continuous group homomorphisms to C\{0} are called 
the characters of G. The group of characters 


G' := {x ■ G —)■ C\{0} : X is a character of G} 


(54) 


is called the dual group of G, because of the duality relation (G')' = G. 

As will become clear after having discussed some examples, the dual group G' can be interpreted 
physically as a “wavevector space” which parametrizes the radiation that solves the design equations. 


just as the wavevectors ko in (43) parameterize plane waves. 


Finally we remark that the design equations as stated in Dehnition 7.1 are intrinsically abelian, and 


should really only be used for abelian G. Namely, (54) shows that the action of G on simultaneous 


eigenfunctions must be abelian, {gig 2 )Eo = ( 5 ' 25 'i)Eo, since the right hand side in (53) is independent 
of the order of the gi. 

For non-abelian but compact subgroups G of E(3), a generalization of the design equations 
which can yield radiation families on which G acts in a non-abelian way has been worked out by 
one of us, and will be presented elsewhere |,Tnl5] . This may be of interest to analyze structures 
such as buckyballs and icosahedral viruses |CK62] . which are generated by non-abelian discrete 
symmetries. In this case the right notion of characters of G is not given by (54), and has the 
structure of a hypergroup |Lal5j instead of a group. 


8 Plane waves as solution to the design equations 

After having formalized our design criterion for structure-adapted radiation into a set of equations 
(Dehnition |7.1[ ), we can state our insight from Section that plane waves are right for crystals as 
a mathematical theorem. 
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Theorem 8.1 (Plane waves are right for crystals). Let G be the translation group T (see 


(40 )j. Then the solutions to the design eguations (Def. 7.1) are precisely the plane waves 

Eo(x) = ko G n G C^, ko ■ n = 0. 


( 66 ) 


Proof By the character lemma, the function y : G —)• C in (50) must be a character of T. The 
characters are well known to be given by 


yk(a) = k G 


Fix ko G and consider the character y_ko- The first eq. ((i) in Def. 7.1) says that 

Eo(- - a) = for all a G Ml 

Evaluation at x = a gives Eo(0) = e“**^°'^Eo(x), that is to say Eq is of form with n = Eo(0). 


The first of the Maxwell equations in Def. 7.1 holds automatically. The second one holds if and 
only if ko ■ Eo(0) = 0. This completes the proof. 

The proof says that hxing a wavevector ko corresponds precisely to hxing a character y : G —)■ 


m 


the symmetry condition (50). For each hxed character, the solutions of the design equations form 


a complex vector space parametrized by {n G : ko ■ n = 0}. This vector space has dimension 2, 
except in the special case y = 1, where the dimension is 3. 

Let us also look at what happens in the closely related case when G is a two-dimensional translation 
group, 

= {(I|a) : a G M^, a ■ e = 0}. (56) 

Here e is a given unit vector in M^. This group naturally arises as the continuous extension of 
the discrete generating symmetry group Go of a 2D crystalline sheet, an important example being 


graphene. For the group (40), we claim that the solutions to the design equations (Def. 7.1) are 
the plane-wave pairs 


Eo(x) = + n_e 


2(1—2e(S>e)ko-af 


n± G C^, ko ■ n+ = 0, (I — 2e (g) e)ko ■ n_ = 0, 


(57) 


with the corresponding frequency in the design equations given by a; = c|ko|. Note that the two 
wavevectors ko and (I — 2e (g) e)ko are the mirror images of each other with respect to the symmetry 
plane k ■ e = 0. 


We sketch the argument leading to (57). We may assume e = (0, 0, 1). The eigenfunction prop¬ 
erty under the symmetry gives Eo(xi, 0:2, (Ts) = 0 , xs) for some planar wavevector 

(fci, k2) G M^. The fact that Eo must be a bounded solution to Maxwell’s equations then gives that 
the dependence on X3 is also of complex exponential form, Eo(0, 0, x^) = -f The 

remaining assertions are then straightforward. 


9 Twisted waves 


We now look at the case when G is the helical group Tie, (41). This group naturally arises as the 
continuous extension of any discrete generating symmetry group of a helix or nanotube structure. 
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We parametrize the group by the rotation angle 6 about the helical axis e, 


( cos 6 — sin 6 0\ 
sin 6* cos^ 0 j 

0 0 l) 


(58) 


It will be convenient to work in cylindrical coordinates with respect to the helical axis, that is to 
say 


'O' 


'Xi 


r cos ip 


e = ea = 0 , x = 0:2 = r sin(y9 , r e [0, cx)), p G [0, 27r), 2 ; G 


(59) 




In cylindrical coordinates, the action (39) of the helical group on vector fields assumes the following 
simple form: 

(Re|re)E) (r, p, z) = ReE(r, p-O.z-r). (60) 


See Figure 3. Here and below, E is a function from polar coordinate space to the cartesian space 


e 



Figure 3: Action of the helical group on vector fields and design of twisted waves. Start from a given 
vector field E. The transformed field E' (right hand side of eq. (60)) under the action of a typical 
element of the helical group is obtained as follows: (i) translate the base points along the helical 
axis e by some amount r while leaving the held direction unchanged; (ii) rotate the base points by 
some angle 0; and (hi) rotate the held direction by the same angle. The design equations require 
that the transformed held only dihers from the original held by a phase factor, regardless of how r 
and 6 are chosen. Time-harmonic electric helds with this property are named twisted waves., and 
can be proven to have the mathematical form (61) (see Theorem 


9.1). 


that is to say Ei, E 2 , E^ are the cartesian held components of E and the action of the group 
on the direction of the held vectors is the usual action of the 3x3 matrix R^ on vectors. 


Theorem 9.1 (Twisted waves). Let G be the helical group "He with axis e. Then the solutions 


to the design eguations (Def. 7.1) are precisely 


Eo(r, p, z) = R^ N(n) 


'Ja+lM 

Ja-lM 

JaM 


(61) 


with (a,/3, 7 ) G Z X M X (0, cx)) (“parameter vector”), n G (“polarization vector”), and 


( 0 , 7 ,/?) -11 = 0 . 


(62) 
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Here N(n) is a certain 3x 3 matrix ( “polarization tensor”) which depends linearly on the polarization 
vector n, 

ni-in2 g 

N(n) = ^ 0 I , (63) 



0 




the Ja are Bessel functions, {r,(p,z) are cylindrical coordinates (59) with respect to the helical axis, 
and Eq is the cartesian field vector. The associated frequency in the design equations is given by 
u = c|( 0 , 7 ,/?)|. 

We call the electric fields (61)-(63) twisted waves. Figurein the Introduction shows the twisted 
wave with parameter vector (a, fi, 7 ) = (5, 3,1) and polarization vector n = (1, 0, 0). 

As already remarked in the Introduction, twisted waves consist of four factors: a scalar plane 
wave on the cylinder; a rotation matrix which rotates the field direction along with the base point; 
a somewhat mysterious polarization tensor; and a vector of three Bessel functions of neighbonring 
order. 

A snbtle aspect is that each twisted wave possesses two analoga of the plane-wave wavevector 


ko in (55): first, the parameter vector {a,l3,'y) which parametrizes the radiation and consists of 
an angular, an axial and a radial wavennmber; and second, the cartesian vector ( 0 , 7 , fi) which 
determines the allowed plane of polarization vectors as well as the freqnency u of the wave (or, 
eqnivalently, its wavelength, see eq. (10)). The latter vector has a simple bnt non-obvious physical 


meaning which will emerge in Section 14 


The parameters a and fi can be interpreted as eigenvalnes of angnlar momentum respectively 


momentum. Namely, it is easily checked that the twisted wave (61), (63) is an exact solntion to the 
eigenvalne eqnations 


J^Eo = a Eo, P^Eo = fi Eq, where J„ 


1 d 





P. 


i dz 


(64) 


Here is the well known qnantnm mechanical momentum operator in axial direction, and Jz is 
the correctly defined angnlar momentum operator on vector fields with respect to the helical axis 
(whose cartesian form (65) can be fonnd in |CT97j b These operators arise in our context of classical 
electrodynamics as the infinitesimal generators of the action (39) of the rotational and translational 
snbgronp of the helical gronp (58) on vector fields: in cartesian coordinates. 


d 

d 

dr 


R0E(Rg ^x) 
E(x — re) 


e ■ (x A V) -I- eAjE(x) = 
(e . V)E(x) = i(rE)(x). 

I 


-(^.E)(x), 

i 


(65) 

( 66 ) 


We remark that the eigenvalne eqnations (64) are eqnivalent to the design equation (50) with 
character x( 6 ',r) = e“4“®+/3'r)_ 


Proof We split the proof into three parts, dealing in tnrn with the three conditions in Def. 
the symmetry condition, the Helmholtz equation, and the divergence condition. 


7.1 


Step 1: Symmetry condition. 

character of the helical gronp "He- 


By the character lemma, the fnnction y in (50) must be a 
The parametrization (58) shows that the helical group Tie is 
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isomorphic to x M, where = [0, 27r) with the usual addition of angles modulo 27r. The 
characters x : *5'^ x M —)■ C\{0} of x M are well known to be 

= («,/?) eZxM. 


In particular, the dual group (54) is given by 




Z X M. 


Fix (a,/9) G Z X M and consider the character X-(a,p) 

((i) in Def. [^) says that 

R,Eo(r, ip-9,z-T) = e*(-“^-^")Eo(r, v^, z). 


X M —)■ C\{0}. The first design equation 


Evaluation ai ip = 9, z = t gives 


Eo(r,v?,^) = e*(“^+^^)R^Eo(r,0,0). 


Step 2: Helmholtz equation. Substituting the ansatz 


(67) 


into the Helmholtz equation (Erst 
of the Maxwell equations in Def. 7.1) gives an ordinary differential equation for Eo(r, 0,0). But 
unlike in the case of the 2D translation group 7^±, due to the presence of the rotation matrix R<^ 
this ODE will not decouple into independent ODE’s for the components. 

The idea to overcome this difficulty is to simultaneously diagonalize the matrices R<^, p G [0,27r), 
with a unitary transformation. Simultaneous diagonalization is possible because the group S'0(2) 
of these matrices is abelian. We have 






^-lif 


U ^ for all p G [0, 27 r), with U = 


V 2 ^/2 ^ 

OOF 


( 68 ) 


Let E(r, 93 , 2 ;) := U ^'Eo{r,p, z). By ( |67| ) and ( | 6 ^ , we have 

p<p 


E(r, p, z) = 


|E(r,0,0). 


(69) 


Since the Helmholtz equation is invariant under the transformation Eq E = f/ ^Eq, we can 
substitute E into this eq. and obtain, using that the Laplacian in cylindrical coordinates is 

^2 Id d^ 

dr'i r dr dp"^ dz^ ’ 

the following ODE for the components: 


^ _ (g + 


+ ((P "-0 = 7 = 1.2,3, 


(70) 


where cxi = 1, (72 = —1, o's = 0. That is to say, the radial functions Ej{r,0,0) are solutions to 
Bessel’s equation. Boundedness of Ej implies that we must have {oj/cY — /?^ > 0. Moreover, for 
integer values of a there is only a one-dimensional space of bounded solutions, given by 


Ej{r, 0, 0) = cj J„+^^.( 7 r), 7 



c,gC. 


(71) 
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By applying U we obtain the original field Eo(r, 0, 0), 

/ J«+i(7r)\ 

Eo(r,0,0) = U diag(c) Ja-iiir) , c G 

\ ) 

with 7 as in ( [7l| . For reasons that are not apparent at this point but will emerge later, it is fruitful 
to parametrize the solution space not by c G but by 


n := Uc 


\ V2 

\ C3 


eC®. 


(72) 


It follows that 

/ J„+i(7r)\ 

Eo(r,0,0) = N(n) J„_i(7r) , n G C^, N(n) := Udmg{U-^n). (73) 

V JaM ) 


Here and below, diag(c) denotes the diagonal matrix whose diagonal entries are given by the 
components of the vector c. Using the explicit form of U (see eq. (68)) it is easy to check that the 
3x3 matrix N(n) introduced above is given by the expression in the theorem. Substitution into 
(67) shows that Eq has the form (61), except that n G is still arbitrary. 


Step 3: Divergence condition. It remains to analyze the second Maxwell equation, divEo = 
0. This is not straightforward, due to the fact that the field components are cartesian but the 
coordinates are cylindrical. We begin by expressing the field components in cylindrical coordinates 
also. 

We denote the unit vectors in the direction of r, if and z by e^, e<^, e^, that is to say 



(74) 


The cylindrical field components are ■ Eq, • Eq, ■ Eq. Since the cylindrical 



unit vectors are the image of the cartesian unit vectors under the rotation R<o, i.e. 




e^f, = R;pe 2 , Bz = Ri^Cs, the rotation matrix cancels from the cylindrical field components: 

E, = E(“‘^+^")Eo(r, 0, 0) ■ ei, E^ = E(“^+^")Eo(r, 0, 0) ■ ea, E, = E(“^+^")Eo(r, 0, 0) ■ eg. 

Moreover, in terms of c = the polarization tensor is 



It follows that 

'Er 

E^ I (r, (p, z) = E(“‘^+^^)N(n) 
Ez 


N(n) = Mfe 


'J„+i(7r) 

Ja_i(7r) 


(75) 


_ gi(o(,5+/32) 


'^(CiJ„+i(7r) - C2Ja-iMy 
^{CiJa+lM + C2Ja-lM) 
CsJaM 


(76) 
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(77) 


Applying the formula for the divergence of a vector held v in cylindrical coordinates, 

dvr I I dv^ dv^ 

div {Vrer + = — + -Vr ^ - 

(jr r r ocp oz 


gives 


divE„=e«“^+'’'> (ci 7 J'+,- C 27 J'_i+ 7^+1- 7^-1) + 7 


zaci iac2 \ , -o T 

- Ja+1~\ - Ja-1 ) +'l-PCsJc 


Here we have dropped the argument yr of the Bessel functions. By Bessel’s identities 






the terms involving Ja+i and Ja-i cancel and we obtain the hnal result 

Cl + C 2 


div Eo = ^ 7 


72 


+ /?C 3 J„(7I-) = i ((0,7, P) ■ n) 4 ( 71 -) 


(78) 


(79) 


In the last equation we have re-expressed the vector c in terms of n. In particular, we see that the 


held is divergence-free if and only if eq. (62) holds. This completes the proof of Theorem 9.1 


Next we present an interesting algebraic property of the somewhat mysterious polarization tensor 
which emerged from the above proof. 


Lemma 9.1 (Intertwining lemma). The polarization tensor ^{n) introduced in (73) satisfies 




R^N(n) = N(n) 


^-vp 


, for all n G C^. 


(80) 


Eq. (80) means that N(n) “intertwines” the standard representation and the diagonal representa¬ 
tion of the rotational subgroup SO{2) of the helical group Tie on C^. 


Proof By the diagonal representation (68) of R(p and the fact that diagonal matrices commute, we 
have, abbreviating the diagonal matrix in the lemma by 

R^N(n) = UD^U-^N{n) = U D^diag{U-^n) = 7/diag(t/-^n)T)^ = N(n) D^. 

The lemma immediately implies an equivalent and sometimes useful representation of the twisted 
wave (61) which we note for future reference: 

Ja+lM 

Eo(r,(p,z) = I Ja-iM 

JaM 



(81) 


Finally we compute the magnetic field associated to a twisted wave. We use the cylindrical com¬ 


ponents (|76|) of the wave and apply the formula for the curl of a vector held v in cylindrical 


coordinates: 


curl {vrGr + + V^efij = 


-e,. 


-7 


JT JL JT 

dr dp dz 

Vr r Vz 
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After some calculation, we find using (78) that 


'j„+i(7r)\ /O' 

curlEo = x n) | , ko ;= 7 

JaM j \Pj 


(82) 


That is to say, to obtain the curl one just has to replace the vector n inside the polarization tensor 
by zko X n. 

The magnetic held associated with the twisted wave (61) can be immediately read off from ([^ 
and (82): 


'J„+i(7r 

Bo{r,(p,z) = R^N(;^ko x n) ( J„_i( 7 r) ) , with ko as in (|82|). 

Mir 


(83) 


Hence for twisted waves, with the “right” dehnition of wavevector the map from Eq to Bq is precisely 
the same map on the polarization vector n as for plane waves (|^. A deeper understanding of this 
fact will be achieved in Section I 


10 The reciprocal lattice of a helical structure 

Crystals are orbits of a hnite set of atoms under a discrete group of translations (a lattice). Helical 
structures are orbits of a hnite set of atoms under a discrete subgroup of the helical group. For 
the latter structures, in the context of X-ray hber dihraction a notion of “reciprocal helical lattice” 
has been introduced by King, Crick, and Wyckham |KCW58j . by periodically extending the helical 
subgroup to a Bravais lattice in and applying the concept of reciprocal Bravais lattice. We show 
here that this notion of reciprocal helical lattice has an intrinsic group-theoretic meaning which 
parallels, rather than needs to rely on, that of the reciprocal lattice in the crystal case. This group- 
theoretic meaning will be very helpful in understanding and interpreting the diffraction patterns of 
helical structures subjected to twisted waves. 


Recall from (42) that the atomic positions in a crystal structure. 


S = -|- a : a G £, = 1,.., M}, 

are the orbit of a hnite set of positions Xq^^ under a Bravais lattice C (i.e., a discrete subgroup 

of the translation group of form C = Al? for some invertible 3x3 matrix A). Analogously, the 
atomic positions in a helical structure, 

5 = {(7x(") : geHo,u = l,..,M}, (84) 

are the orbit of a hnit set of points Xg^\ under a discrete helical group, by which we 


mean a 


discrete subgroup of the helical group "He with axis e (see (41)) of form 


^0 ^ {OoK ■ * £ j ^ C = (R(>o|roe), 50 ^ (R2Vr<|0), ^0 e [0,27r), To € M, To # 0, n e N. 

( 86 ) 

In case n = 1, 5^0 is the identity and 5 is a helix. In case n > 1, 5 is the union of n helices, 
basic examples being “zigzag” or “armchair” carbon nanotubes. The parameters Tq and 6*0 of the 
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generating screw displacement Hq encode the pitch p and the number q of subunits per turn of the 
helices. The pitch is dehned as the axial displacement for one full rotation, 


2ti 


( 86 ) 


and the number of subunits (i.e. rotated and translated copies of the set {xg^\ Xq^^}) per turn is 

27r 


q = 




(87) 


In physical and biological examples, the latter number is typically a rational number but not an 
integer. 

Before giving explicit examples, we bring formula (85) into a form more similar to a Bravais 
lattice, by writing the group elements in the form (R^lre). We have 

. . 27r 

5'o^o = (R-el're) with 6 = —i + 6*ojmod27r, r = tqJ. 

(Here 6*mod27r denotes the unique angle in [0, 27r) which differs from 9 by an integer multiple of 
27r.) The above equations for 9 and r can be written more compactly as 


T 


= I 


2-k 


0 


+ J 


mod 


27r 

0 


^ 0c 

n ^ 

0 tq) \j 


mod 


27r 

0 


Denoting this parameter set by T-Lq, i.e. 

Tio = AZn X Z mod 

it follows that 


27r\ . , . 9t 


with A = 


0 To r 


Ho = {(Re|re) : 


G Ho}, 9o G [0, 27r), To G M, To 7 ^ 0, n G N. 


( 88 ) 


(89) 


(90) 


(Here, nmod(27r,0) denotes the unique vector in [0, 27r) x M which differs from v by an integer 
multiple of (27r,0).) See Figure]^ 


Example 1: Carbon nanotubes. All nanotubes are of form (84), (90). To give a specihc 
example, single-walled (6,5) Carbon nanotubes with axis e = (0, 0,1) correspond to 


M 

= 2 (number of atoms per unit cell) 

n 

= 1 (single helix) 


Oo 

182 


To 

= —= i with i = 

2yM 

1A3A° (C-C bond length) 

x(i) 

= (r, 0,0), where r 

= i (nanotube radius) 

ZTT 

X(2) 

= (R„,/3lie)x«. 
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Example 2: Tobacco mosaic virus (TMV). This is a basic example of a filamentous virus, built 
from a single protein. The protein molecules are arranged in a low-pitch helix, with 161 proteins 


per turn and with adjacent turns in contact. The parameter values in (84), (90) are 


M 

n 

Oo 

To 


= 1284 (number of atoms in the protein) 
= 1 (single helix) 

27r 3 

= —Y = — 27r 

16| 49 

= 1.402 AT 


The values for 9q and tq, taken from |GZ11] . correspond to the low-calcium state. 


The groups in (90) are not some ad hoc ansatz. It can be shown that they are the lattice subgroups 


of the helical group (41), whereas the Bravais lattices are the lattice subgroups of the translation 
group (40)]^ 


We now introduce a purely group-theoretical notion of reciprocal lattice. We first state this 
notion in abstract mathematical language, then show how it reduces to familiar concepts in the 
crystalline and helical case. 

Definition 10.1 (Reciprocal lattice group) Let Hq he a lattice subgroup of the helical group 
(see (90 )j, or more generally any lattice subgroup of a locally compact abelian group H. Recall from 
(54) the dual group of H, H' = {x '■ R —)■ C\{0} : y is a character of H}. The reciprocal lattice 
group of Ho with respect to H is the set of those characters of H which are egual to 1 on Ho, i.e. 


H'q = [x & H' : x{h) = 1 for all h G Ho}. 


Thus mathematically, the reciprocal lattice group H( is a subgroup of the dual group H'. Physically, 
the dual group consists of certain scalar waves parametrized by wavevectors or wavenumbers (see 
the examples below), and the reciprocal lattice group consists of resonant waves, parametrized by 
a subset of wavevector or wavenumber space which we call the reciprocal lattice. 

In group theory, H( is a well known object, called the annihilator of Hq (see e.g. |HR63] j or the 
orthogonal group of Ho (see e.g. |RS00j : this terminology views the equation x{h) = 1 as analogous 
to the vanishing of a proper inner product between elements y and h of some vector space) with 
respect to H. 

Example 1: Reciprocal lattice of crystals 

a) Mathematical description. Let H be the translation group T = {(I|a) : a G M'^}, and let C be 
a Bravais lattice, i.e. C = Al/, A an invertible d x d matrix. The characters of H as functions of 
a G are known to be the scalar plane waves 


Xk(a) = e 


ika 


k G 


tid 


An group theory, a subgroup is called a lattice subgroup if it is discrete, i.e. has no accumulation points, and 
there exists a set of finite volume (Haar measure) whose orbit under the subgroup gives the whole group. To derive 


the representation (85), one shows and uses that each such subgroup must be the image of a lattice subgroup of the 
two-dimensional translation group under the map —>■ He defined by h{ 9 ,T) — (Remod27rl'’'®)- 
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Hence dual group H' and reciprocal lattice group Hq are 

H' = {Xk : k G M"}, 

Hq = {xk : k G k ■ a G 27iZ for all a G H Z'^}. 

b) Physical description. The characters Xk ^ H' are parametrized by the wavevectors k G The 
wavevector space which parametrizes the group H' is 

PL' = 

and the set of wavevectors which parametrizes the reciprocal lattice H'q is 

PL'q = {k G : k ■ a G 27rZ for all a G 


recovering the standard dehnition (34) of the reciprocal lattice. 

Example 2: Reciprocal lattice of helical structures. 

a) Mathematical description. Let H be the helical group {(Rglre) : 6 G [0,27r), r G M} (see (58)), 
and let Hq be a lattice subgroup, i.e. a group of form ( pOj ). The characters of R as a function of 6 
and r are known to be the cylindrical waves 


Xa,/3(^,r) = e*(“^+^"), (a,/?)GZxM. 
Hence dual group H' and reciprocal lattice group Hq are 

= {Xo,/! : (“./3) 6 Z X R}, 


(91) 


(92) 


K = {X.,,11 ■■ (a,;3) 6 Z X R, (^“j • (^^j € 27rZ for all (^^) E AZ„ x Zmod (^J") }. (93) 

b) Physical description. The characters (91) are parametrized by an angular wavenumber a G Z 
and an axial wavenumber /3 G M, yielding the following parameter space for the dual group H'\ 


H' = {{a, 13) : oGZ,/?gM} = Zx 


The reciprocal lattice is 


U'q = {{a, (3) G Z X M : 


G 27rZ for all 


r 


G x4Z„ X Zmod 


27r 

0 


(94) 


}. (95) 


Using the explicit form (88) of Ho, it follows that (a, (3) belongs to the reciprocal lattice if and only 
if 

27r 

a — G 27rZ, aOo + (3tq G 27rZ. (96) 


n 


The solutions (a,/?) to this condition are easily computed and one obtains, analogously to ( piSj ), 
PL'n = 2'kA~'^Z X Z = 


nOp 2n ) Z X Z 

TO TO 


ndp j + J f ^ ) ■ hJ ^ 

TO / V TO , 


(97) 


Formula (97) recovers the reciprocal helical lattice associated with a helical structure generated by 
the group (90) as introduced in |KCW58] . and reveals its group-theoretic meaning as a wavenumber 
space which parametrizes the group (93). See Figure]^ 
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Figure 4: A discrete helical group (left, eq. (90)) and the associated reciprocal helical lattice (right, 


eqs. (93), (97)). Here 6 is the rotation angle about the helical axis and r the displacement along 
the axis, and the reciprocal parameters a and f3 are an angular respectively axial wavenumber. The 
slope of the reciprocal basis vector pointing to the right is the inverse pitch multiplied by 27r, by 


eq. (86). The discrete helical group is a subset of the full helical group 'He=*S'^x M (eq. (58), shaded 
vertical strip), whereas the reciprocal helical lattice is a subset of the dual group Z x M (eq. 
shaded vertical lines). 


11 Diffraction of twisted waves: twisted Von Lane condi¬ 
tion 


We now calculate the outgoing radiation when twisted waves are scattered off helical structures 
of inhnite length. The reciprocal helical lattice will naturally appear, due to its intrinsic group- 
theoretic meaning derived in the last section which leads to an associated Poisson summation 
formula on the helical group. 


Let Eo(x;a,/3,7) be a twisted wave (61) with parameters (a,/?,7) G Z x M x (0, cx)). By eq. 
(^-(ig), the outgoing field scattered from a structure with electron density p is, assuming Xc = 0, 




32 (k(x)-x—cjf) 


Ik'(x) 


Ik'(x) 


fEo(k'(x)) 


with outgoing wavevector k' and structure factor Ieq given by 


k'(x) = fEo(k'(x)) = 

C X 




(98) 


(99) 


The difficulty is to evaluate the structure factor and unearth the fact that it has a sharp peak 


structure, which is far from obvious from eq. (99). 


We will proceed in 5 steps: (1) transform to cylindrical coordinates; (2) eliminate the ensuing 


phase nonlinearity of the plane-wave factor inside the integral (99) by expanding this plane wave 


into cylindrical waves, which mathematically corresponds to a Fourier series expansion in the angle 
variable (this step is not necessary in case of an axial detector); (3) treat the integration over (p' G 
and z' G M jointly rather than separately, because these variables are “intertwined” in the electron 
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density of a helical structure, and use Fourier calculus on the group x M to essentially reduce 
the integral to the Fourier transform of an infinite sum of delta functions on x M; (4) apply the 
Poisson summation formula on the abelian group x M, thereby obtaining a twisted analogue of 
the Von Lane condition in case of an axial detector; (5) eliminate the expansion of Step 2 in case 
of a non-axial detector. 

Step 1 (Cylindrical coordinates). Let {r,(p,z) and {r',Lp',z') be cylindrical coordinates for x re¬ 
spectively y, i.e. X is given by (59) and yi = r' cos ip', 1/2 = r' sm.ip', = z'. Next, let (i?, <F, Z) be 

the cylindrical components of the outgoing wavevector held k'(r, ip, z), i.e. k' = Ror + -|- Ze^, 

( 100 ) 


with the cylindrical unit vectors of eq. (74). It follows that 


R{r,ip,z) = - T===, ^{r,p>,z) = 0, Z{r,p>,z) = 


c \/r 2 -f 2 ) 2 ' 


c 


With the help of the rotations and (see (58)), we have 


k'(r,(p,z) = R,^ I 0 I , 
Z. 


y = R^/ 0 . 


Consequently the phase of the plane-wave factor in the integral (99) is 

k' ■ y = Rr' cos{ip — ip') + Z z'. 


Substitution into the representation (81) of twisted waves in which the rotation matrix has been 
diagonalized gives 




fEo(k'(r,v?,^)) = N(n) 


' ir! 


/o 9SlxI 


p(r', (p', V) e 


o-w 


'j„+i(7r') 
J„_i(7r') 
J„(7r') 




( 101 ) 


where here and below, by J^^^^^dip'dz' we mean dip'dz'. 

Step 2 (Eliminate the phase nonlinearity). We expand the plane wave factor in (101) into cylindrial 
waves, corresponding to a Fourier series expansion in the angle. The Fourier expansion of the 27i- 
periodic function is the Jacobi-Anger expansion 


—iA cos 9 


Y, {-irMA)e’-> = 5 ^ Z{A)> 


,^1^(0—7r/2) 


( 102 ) 


where the J,, are Bessel functions. It follows that 


/•C.XJ 

fEo(k'(r,v^,^)) =N(n) / J.(Rr') 


+{I3-Z)z'] 


'51x1 


p(r', ip', z' 


'Ja+i{ir') 


( 103 ) 


g«[(o 1 vW+iP z)z'] dip'dz' j jQ_i( 7 r') r'dr'. 


p[{a-v)ip'+{P-Z)z'] 


Jai'yr') 
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In the special case of an axial detector, i.e. {r,ip,z) = (0,0, z), the cylindrical components of the 
ontgoing wavevector k'(r, y?, z) are 


u 


R = 0, d) = 0, Z = — sign( 2 ;). 


( 104 ) 


Thns the phase nonlinearity is absent and we obtain directly from (101) that 


fEo(k'(0,0,z)) = N(n) 


0 '-Jgixi 


p(r',(p',z') 


Ja+liir') 

^i[{a-i)ip'+{i3-z)z'] I I jQ_i( 7 r') I r'dr'. 


J[aip'+{P-Z)z'] 


Jahr') 


Note that in the case (104), the sum (103) indeed reduces to (105) as it should, because of the 
property of Bessel functions that 

Ju{0) = 1 for z/ = 0 and 0 otherwise. (105) 

Step 3 (Exploit helical symmetry and Fourier calculus on x M). Consider now a helical structure, 
i.e. a structure generated by any discrete helical group Hq (see (90)). The electron density p will 
be ifo-periodic, that is to say 


p{h y) = p(y) for all y G and all h G Hq, 


(106) 


and rapidly decaying in the direction perpendicular to the helical axis. Typical examples of Hq- 
periodic densities are depicted in Figures 6 and 7. By choosing a suitable partition of unity, such a 
p can be written as a sum of rotated and translated copies of a localized, rapidly decaying function 


p(y) = 

heHo 

In cylindrical coodinates, iifo-psriodicity means that 


( 107 ) 


p(r', ip', z') = p{r', if' — 9 mod 27r, z' — r) for all 


G Ro, Rq = AZn X Zmod 


27r 


and the representation (107) means that 


p{r', ip', z') = — Oi mod 27r, z' — 02 ). 

aS'Ho 


( 108 ) 


The function can for instance taken to be the restriction of p to the unit cell U of the structure. 


= pu = 


p in 7/ 


U = {(r', ip', z') : {ip', z') G A [0, 1]N r' G (0, cx))}. 


(109) 


0 outside 7/ ’ 

but other constructions with a smooth 'ifj make sense too. 

The decomposition (107) of p can be fruitfully re-written as a convolution of 'ip with an inhnite 
sum of delta functions, 

p — p) ^T-Lo — '^a, 











where the convolution on x M is dehned as 


if = / f{r,ip-ip'mod27r,z-z')g{r,ip',z')dip'dz'. 


'Sixl 


( 110 ) 


In the sequel we drop the subscript from the convolution sign. We now use Fourier calculus on 
X M. The Fourier transform with respect to the angle and axis variables is a function on 
the dual group (5^ x M)' = Z x M, dehned as 




2ti 


( 111 ) 


'S^xl 


The square bracket in (103) has the form of such a Fourier transform on x M: 

p^,^^(V^*<5«J(r',i.-(a + l),Z-/?) 0 0\ 

= 27r 0 (a-1),Z-/?) 0 . 

Vo 0 Jl,Ji;*6u,)ir',u-a,Z-fi)J 

Step 4 (Poisson summation formula on x M). We now use the (trivial) convolution rule on 
X M, 

jr = 271 / ■ -F, g, (112) 

and the following nontrivial result from Fourier analysis on abelian groups which makes the recip¬ 
rocal helical lattice appear: 


Lemma 11.1 (Poisson summation formula on x 


Let Hq be any discrete helical group 


(see eq. (90 )), and let H'q be the reciprocal lattice. The corresponding parametrizations TLo G x M, 

2% . 


(89), and TL'q C Z x M, (97), satisfy 




’n'a 


(113) 


This identity is a special case of the general Poisson formula on locally compact abelian groups 
going back to A. Weil |We64] . see e.g. |RS00j ^ A more elementary derivation of (113) is to hrst 
consider the case 6^ = 0, where the result follows by combining the usual Poisson formula on M with 
the following, elementary to check, Poisson formula on S^: if F = {— : j = 0, ..,n — 1}, then 


n—1 


SM = T e [0,27r), Sciiy) = ^ tz e z. 


j=0 


aGnlj 


in 


Here f{u) denotes the Fourier coefficient (27r)“^ f (t) Note that the delta functions 

the left sum are Dirac deltas, whereas the delta functions in the right sum are Kronecker deltas. 
The general result (113) now follows from a suitable change of variables. 


Eqs. (112), (113) yield 


* *^^ 0 ) I det A| 


( 114 ) 


®In this context, the formula is stated and derived up to an overall multiplicative constant. 
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and therefore 


fEo(k'(r,(^,z)) = 


det A\ 


“ - {(A-l),Z -/3y 

N(n) Z) Sn;,(i' - (cH-1), Z-ff)], (115) 

u=-oo \ Sn'gii^ - a, Z -/3) 


with Z) = I Ju{Rr') ( 0 - (a - 1), Z - (3)Ja-i{'^r') 0 | r'dr'. 

‘0 0 


For an axial detector (104), the sum over u reduces to the contribution from the single term i/ = 0, 
thanks to (105). Moreover for u = 0 the hrst Bessel factor, Jy{Rr'), in the matrix-valued integral 
in (115) equals 1. Hence the matrix components of the integral involve just one remaining Bessel 
factor, and can thus be interpreted as a Hankel transform. Recall that for any a G Z, the Hankel 
transform of order a maps scalar functions of a radial variable belonging to the interval (0, oo) to 
scalar functions on (0, oo), and is dehned as 


{Haf){l)= !(r')Ja('ir')r'dr' (7 > 0). 


(116) 


It follows that 


fE,(k'(0,0,i)) = 


(2i) 


det HI 


^5^i(-(cH-l),Z-/3)^ 
N(n)Zl„,^,,(Z) I 5H'S-{oM),Z-y) 


6n'y-a,Z-y) 

+ - y) 0 0^ 

with = I 0 - y) 0 

,0 0 


(117) 


(118) 


For the purpose of calculating the outgoing intensity, it is useful to write the delta function 5-^^ as 
a sum over reciprocal lattice vectors, 5-^^ = Yha&u'Q ^kat 

5a'(-(« ± 1), ^ - /3) = (-a, Z-y). 

This yields the following alternative expression for the structure factor: 


fEo(k'(0,0,z)) = 


det H| 


a'e-H^ 


1 

X 



V / 


{-a,Z-y). 


(119) 


Note that the delta functions in (119) are centered on shifted copies of the reciprocal lattice. Also, 
we claim that the Fourier-Hankel transform of fj which appears in (119) equals that of the electron 
density (109) in the unit cell, i.e. 

^ ^0- (120) 
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This is because, for any ^|J satisfying (107), including ^|J = pu, the Fourier transform a') 

is independent of 'ijj when a' is a reciprocal lattice vector; indeed 




'51x1 


e ''ip{r',ip', z')dip'dz' 


m 


e **^*'(^0 — Oi mod27r, z'— 02 ) 

e (^', z) dif'dz' for all a' G "Hg. 


>u 


Here we have used that the factor in the middle line equals 1, and have employed (108). 


The outgoing electric field can now be read off immediately from (98), (99), (119), (120). Note 


in particular that the third held component vanishes, since the projection matrix in (98) annihilates 


the third component of the structure factor when the outgoing wavevector k' points in axial direction 
± 03 . See eq. ( 121 ) below. 


When taking absolute values to obtain the intensity, it may happen that the two remaining 


shifted copies T-Lq + (1, 0) and Hg — (1, 0) of the reciprocal lattice appearing in (119) overlap, which 
would lead to interference. These two copies overlap if and only if (2, 0) is a reciprocal helical lattice 
vector. The following result, which is elementary to check, shows that this does not happen except 
in a degenerate case which we propose to denote flat helical groups. 

Lemma 11.2 (Flat helical groups). The following three statements about a discrete helical group 


(85) are eguivalent: 

(1) The vector (2,0) belongs to the reciprocal helical lattice, eg. (97) 


(2) The parameters of the helical group satisfy Oq = 0 or n, and n = 1 or 2. 

(3) The structure generated by applying the helical group to any single point lies in a plane. 

We summarize our hndings as a theorem. 

Theorem 11.1 (Twisted Von Laue condition). Consider a helical structure with electron den¬ 
sity p : —)■ M, assumed to be smooth, HQ-periodic with respect to some discrete helical group Hq 


(see eg. (85)j, and rapidly decaying in the direction perpendicular to the helical axis. Assume that 
the axis is e = Bq, and let the incoming electric field be a twisted wave with same axis and parameter 


vector (a,/3,'y) G Z x M x (0, 00 ) (see (61)-(63)j. Recall that the freguency of this twisted wave is 


oj = c|(0,7,/3)|. Then the diffracted electric field (98) at any point (0,0, z) on the axis is 


Eout{0,0,z,t) = 


i(- 

gH c 

\z\—LOt) 

CO 


(H 

^ei 

z\ 

det A 



/ ni+in2 
/ 2. 

711—2712 

2 

0] 



1 n2—ini 

\ 2 

712+1711 

2 

0 


' 0 ) 

V 0 

0 

0/ 

V 0 

/ 


( 121 ) 


(-«, ysignz - fl). 


Here H^ is the Hankel transform of order a with respect to the radial variable (see (116)j, H 


is the Fourier scries/transform with respect to the angular and axial variables (see (lll)y), pu is 


the restriction of the electron density p to the unit cell (see (109)j, and TLq is the reciprocal helical 


lattice (see (93)J. 


Sl> 
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Moreover, when Ho is not a flat helieal group (see Lemma 11.2), the square root of the outgoing 
intensity is 


J(0,0,z; a,/3,7)) ' = 


det A\ \z\ ^ 


( 122 ) 


\ni + m2| 

~^/r~ 




+ 


\ni - m 2 1 

~^/r~ 




where cq = In particular, constructive interference occurs if and only if the difference 

between the angular/axial part ( 0 , ^ signz) of the outgoing wavevector and the angular/axial param¬ 
eters (a,/3) of the incoming twisted wave belongs to the reciprocal helical lattice shifted left or right 
by precisely one angular wavenumber. Ho ± (J), or equivalently, if and only if 


— a = inlul, Toifl: sign z — fl) = —inOo + 27ij for some integers i, j. 


(123) 


Formulae (121), (122), (123) are the main result of this paper. They say that the signal of a 


helical structure under co-axial twisted X-rays, recorded along the axis, consists of sharp peaks 
with respect to the angular and axial radiation parameters. The peaks are double-peaks with a 
distance of precisely two angular wavenumbers, centered at the reciprocal lattice vectors of the 
helical structure. In particular, at the reciprocal lattice vectors themselves the signal vanishes. The 


structural parameters tq, Oq in (85), or equivalently the pitch and the number of subunits per turn. 


can be immediately read off from the peak locations, as can the order n of any rotational symmetry. 

Moreover the above result makes it in principle possible to determine the electron density p, i.e. 
the detailed atomic structure, from intensity measurements in the far held at a specihc point on 
the axis, provided the scalar phase problem associated with the Fourier-Hankel transform 
can be solved. Note that the values of this transform on the reciprocal helical lattice points which 


appear in (122) completely determine the electron density; see Section 12 for details. 


Next, we report an important invariance property of the axial signal of a helical structure. 
Suppose the structure is translated by an amount Az along the helical axis, and rotated by an 
amount Ap around the axis. (In other words, we apply an arbitrary element of the continuous 
helical group (41) to the structure.) This modihes the original electron density pu to 


pu{r, p, z) = pu{r, p - Ap, z - Az). 


From the dehnition of the angular/axial Fourier transform, (111), it is clear that 


SlxR 


pu') (r, I., Z) = (7.„p„) (r, J., Z). 


(124) 


(125) 


Thus this transform changes just by a phase factor. The Hankel transform with respect to the 
radial variable does not interfere with this phase factor, and so it follows that 




for all 7 . 4. Z. (126) 


In particular, the outgoing intensity (122), which only depends on the absolute value of the ex¬ 
pression (126), is invariant under axial translations and rotations of the structure. This is not 


a fortuitous accident, but stems from the fact that the design equations, Def. (7.1), require the 
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incoming wave to be an eigenfunction of each element of the continuous helical group (41). From 


this one easily sees that the invariance remains true for helical structures of hnite length. 

By comparison, the signal produced by hber diffraction, i.e. by sending plane wave X-rays 
towards a helical structure from a perpendicular direction, is only invariant under axial translation 
but not under axial rotation; the latter is a well known major problem in the interpretation of hber 
diffraction images. 


Step 5 (Arbitrary outgoing direction). Finally, we calculate the outgoing radiation in arbitrary 
direction. Before doing so, we note the simple group-theoretic reason for why one expects the 
radiation to be much more complicated when dropping the restriction to axial outgoing wavevectors. 
Only in this case is the factor 6“**^ appearing in (99) a character of the helical group He, i.e. 
an element of the dual group H' of the helical group, which - as explained in Section - is the 


parameter space for twisted waves. But this is necessary in order that the outgoing wavevector just 


causes a shift of the incoming radiation parameters as in (121)-(122). 


Indeed, the outgoing radiation in non-axial direction is more complicated. We start from the 


expression (115), and focus on the hrst component of the (vector-valued) sum over u. We have, for 
any three functions /, g, h on Z, 


(o^^)) h{cH-l)Su'^{u - {aH),Z - f3) = ^ /(a) + aH) g{a[) 5a'^{Z - f3). 

V 

It follows that 


a'en'o 


/Ja[+a+l{R'^ )Ja+l{'l'^ ) 

fEo(k'(r,(p,2;)) = ——— N(n) ^ (Z -/1)F)„/,„((p) / (J^^^^^/>)(r', a') Ja/+„_i(Hr') Ja-i(7^^') r'dr' 

R \ Ja',URr')Jahr') I 


0 


with = I 0 


0 

0 


0 


0 




Denoting the product = Jy{Rr')'il){r\g)\ z') by / = J^{R-)i/j and using eq. (109), this 

formula for the structure factor can be written more compactly as 


fEo(k'(r,v>,z)) = 


(27r) 


^+a+l{R ■ )Pu)\ 


det A\ 


N(n) 


a’&H'r 


0 


S^xK 

Ha-lR 


si XR 


'Ja'j^+a—1 {R ■ )pu 


\ 


“ SlxR 


'Ja\+a {R ■)pu^ j 


(7,a'). 


(128) 


This is the generalization of formula (119) to an arbitrary outgoing direction. Thus the far held 


observed from any direction still exhibits a sharp peak structure in the axial wavenumber jd. More 


precisely, by (97), resonance occurs when 

tq{Z — (3) = —inOo + 271 j for some integers i,j. 


(129) 


But the signal is spread out over angular wavenumbers a that need not be related to the reciprocal 
lattice. In the limit of the radial component R of the outgoing wavevector tending to zero, the 
Bessel factors Ja'^+a+a{R ■) (<7 = 0,1, —1) in front of the unit cell electron density pu converge to 
the delta functions 6a'^+a+a, reducing (128) back to the axial formula (119). 


(127) 
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If the structure factor is averaged over outgoing wavevecors with hxed angle to the helical axis, 
sharp peaks appear also with respect to the angular wavenumber. More precisely, the average of 
the diagonal matrix Dq/,^(( p) over (p is a diagonal matrix of delta functions, eg. the upper left entry 
is 1 if a'^ + a + 1 = 0 and zero otherwise. As a consequence. 


r- 27 r 


— / fEo(k'(r,(p,z))d(p = 


(2vr) 


27r 


det A\ 


a'en'a 




) (7, a') 

7-(;) 


S.' / 


(-a,Z-l3). 


(130) 


The only difference as compared to the structure factor (119) in axial direction is the zero*^ order 
Bessel factor in front of the unit cell electron density. In particular, by (97), resonance of the 
component of the averaged structure factor occurs when 


— a = in + ak, to{Z — fd) = —i n 9o + 27ij for some integers i, j, (131) 

where = 1 for fc = 1, — 1 for k = 2, and 0 for k = 3. If this averaged signal can be realized 
experimentally, the observation direction (which influences the axial component Z of k') can be 
varied to detect resonance. 


12 Synthesis of electron density 


As shown in the previous section, subjecting a helical structure to twisted X-rays and recording the 
intensity of the scattered radiation in axial direction yields the data set 


{|G'(7,a,/?)P : 7 G (0,oo), {a,/3) G Hq} , 
where G is the Fourier-Hankel transform of the unit cell electron density pu, 

1 


(132) 


G(-ua,l3)= {HaT^^^^Pu)(l,a,l3) = — / e J^etr) pu(r,'p,z)r dr dtpdz (133) 

0 v 'J ZA 

and "Hq Ike reciprocal helical lattice of the structure. Assume that the phases of the Fourier- 
Hankel coefficients G in (132) can be retrieved. Then the electron density in the unit cell is recovered 
by 


Pu{r,g},z) = ^ X] (h^G{- ,a,l3)^{r)e 


,i{a(p+gz) 




1 


J] / J„(7r)G(7,a,/?)7(i7e“”'’+'’'>. 


(134) 


&H' 


To show this, one hrst notes that the Fourier-Hankel transform in (133) is a combination of a Fourier 
series in the angle, a Fourier transform in the axial variable, and a Hankel transform in the radial 
variable, and is thus invertible. One then exploits formula (114) with 'ip = pu, and uses that the 

(135) 


inverse of the Hankel transform (116) is given by 


Hp, f){r)= / f{-f) JaM-fd'y. 
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The overall situation is thus the same as in standard X-ray crystallography: the electron density 
can be reconstructed provided we can solve a scalar phase problem. The only difference is that we 
have to deal with a 


New phase problem. Reconstruct a function pu given the absolute value of its Fourier-Hankel 
transform G defined by (133). 


Numerical investigations which will be reported elsewhere |JJF15] indicate that standard phase 
retrieval algorithms for the Fourier transform can be adapted without difficulty to the Fourier- 
Hankel case. 


13 Simulated diffraction pattern of carbon nanotube and 
tobacco mosaic virus 


Here we present simulated diffraction patterns of helical structures subjected to twisted waves. The 
set-up is as in Figure of the Introduction, that is to say incoming waves, structure and detector 
are axially aligned. As pointed out earlier (see (126)), the signal is invariant under axial translations 
and rotations of the structure. 

We consider two examples, a (6,5)-Carbon nanotube and TMV virus. The parameters are as 
in Example 1 respectively 2 of Section with the TMV atomic positions taken from the Protein 


Data Bank, PDB ID 3J06 |GZllj . We use the electromagnetic model (98), (99), (24) to calculate 
the signal from structures of finite length: 


• C nanotube: 255 unit cells, corresponding to 510 atoms. 

• TMV : 147 proteins, corresponding to 3 helical repeats and 188 748 atoms. 




0 


a 



Figure 5: Twisted X-ray pattern of a (6,5)-carbon nanotube with 510 atoms. The electron density 
of each atom was modelled by a Gaussian. Right: Iso-surface of electron density. Top: Log 
intensity of diffracted radiation as a function of angular wavenumber a and axial wavenumber (3 of 
incoming wave. Bottom: Theoretical peak locations, eq. (122). Incoming twisted X-ray wavelength: 
A = 1.54 A° (Cu Ka line). 
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Figure 6: Twisted X-ray pattern of 3 helical repeats of TMV (147 proteins, 188 748 atoms). For 
simplicity, atomic form factors were ignored and the electron density of each atom was modelled as a 
delta function. Right: Iso-surface of electron density. Top: Log intensity of diffracted radiation as a 
function of angular wavenumber a and axial wavenumber /3 of incoming wave. Bottom: Theoretical 
peak locations, eq. (122). Incoming twisted X-ray wavelength: A = 1.54 A° (Cu line). 


The results in Figures 5 and 6 show that already relatively short pieces of a helical structure 
exhibit the theoretically predicted patterns. In particular, one sees double-peaks centered at the re¬ 
ciprocal helical lattice points which would be difficult to interpret correctly in terms of scalar models 
of the incoming wave and the diffraction pattern. Moreover, very high contrast between peaks and 
background is observed; note that in fiber diffraction the contrast in the direction perpendicular to 
the fiber is limited for theoretical reasons |CCV52j . 


14 Fourier representation of twisted waves 


We now derive the expansion of twisted waves into plane waves, or mathematically: the Fourier 
decomposition of twisted waves. This leads to a simple physical interpretation of the polarization 


vector n and the “cartesian reduced wavevector” ko = (0,7,/3) appearing in (61)-(63), both of 


which have remained somewhat mysterious up to now. It also suggests a possible route to realizing 
approximate twisted waves experimentally. 

Consider, e.g. the axial component of a twisted wave, const Jaijr). This expression 

is a scalar cylindrical harmonic, evaluated at the point with cylindrical coordinates {r,(p,z). The 
cartesian coordinates of this point are Xi = rcos<p, X 2 = rsin:^?, x^ = z. Using the Bessel integral 


r-27r 


J^[A) = — / 
27r Jo 


gives 


1 1 p27V 

_ _ / ^i[a(ip'+ip)-ir sirup' 

27r Jo 2 ti Jo 


(136) 


where we have employed the substitution ip' -t- (p = 4). The phase appearing in the latter integral 
can be interpreted as a cartesian inner product. We have, using the rotation matrices and R$ 
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(see (58)) and the addition formula for the sine, 


7 sin$ 

7 r sin($ — ip) + fiz = \ 7 cos $ 

/S 



= R$ko ■ X, 


(137) 


with ko = ( 0 , 7 ,/ 3 ). Eqs. (136)-(137) allow to express the twisted wave (61) as a function of the 
cartesian coordinate vector x. Starting from the equivalent expression (81) in which the rotation 
matrix has been diagonalized, we obtain 


V[(“+i)‘^+^^lj„+i(7r)' 

Eo(x) = N(n) I J„_i( 7 r) 1 = N(n) — 

e*[“^+^^]j„(7r) ' 


-271 




(138) 


The geometric meaning of this expression becomes clear by using the following identity which one 
obtains by applying the intertwining relation (80) to the vector ( 1 , 1 , 1 ): 

N(n) ' = R$n. (139) 

This gives 



-271 


Eo(x) = ^ 






(140) 


Expression (140) is a Fourier reconstruction integral over a one-dimensional circle in k-space. 


C = {R#ko : $ e [0, 2ti]} = {k G 


a; 


■ = ~’ k-e = ,3}, 


(141) 


where e is the helical axis. Here we have used that |ko| = u/c. See Figure 7. 

Since the circle has radius 7 , the line element (or one-dimensional Hausdorff measure) on C 
is ds = 7 (i$. Dehne the following polarization vectorheld as a function of wavevector which 
corresponds to rotating the polarization along with the base point: 


n : C —)■ C^, n(R$ko) = e*“'*’R$n(ko), n(ko) = n. 


(142) 


In terms of this vectorheld on reciprocal space, eq. ( |140 ) can be written in cartesian reciprocal 
coordinates as 


Enfx) = 


1 f (27r) 

(2vr)3 y^ec 7 


-fi(k)e'‘^-"ds. 


(143) 


This is a Fourier reconstruction integral, and hence the Fourier transform of a twisted wave is 


Eo(k) = —^n(k) ds 


7 


(144) 


Mathematically, this is a singular measure supported on a circle in wavevector space. 

In cylindrical coordinates {R, <h, Z) in reciprocal space, the Fourier transform of a twisted wave 
acquires the simple form 

(27r)2 


Eo(R,$,Z) = - ^) 6{Z - /3). 


7 


(145) 
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Figure 7: Fourier transform of a twisted wave. The circle of wavevectors and associated polarization 
directions in the picture indicates the plane waves whose phase-modulated superposition gives the 
twisted wave with parameters (a, /3, 7 ). The fact that just a one-dimensional family is needed means 
that a good approximation can be achieved by a relatively small number of plane waves. This might 
be of interest for the experimental realization of twisted waves. 


Here the radial delta function is, as usual, normalized with respect to RdR not dR, that is to say 
J^6{R--f)RdR=l. 

Formula ( 140[ ) has an interesting alternative interpretation, namely as a group average instead of a 
k-space integral. Start from the plane wave 

Epz,„e(x) = ne**^“-", (146) 

and note that the phase and the plane wave in (140) satisfy R$ko ■ x = ko ■ Rj^x, R$ne**^“'^ = 
R>i.Epia„e(R$^x). Consequently 


Eo(x) J Xaid) \ deplane 


(x)d/i(^). 


(147) 


where G is the group of rotatons around the cylindrical axis, Xa{g) = is a character of the group, 
and dpi = ^d^ is the Haar measure on the group, normalized according to the usual convention 
for compact groups that = 1. This shows that a twisted wave is an integral of the image of 

a plane wave under the group of rotations about a fixed axis against a character. The polarization 
vector n and the “cartesian reduced wavevector” ( 0 , 7 ,/d) are just the polarization vector and the 
wavevector of this plane wave; the angular wavenumber a comes from the character. 

Eq. (147) is an example of a Wigner projection, hrst introduced in the context of quantum 
systems in |Wi31] . Such a projection maps solutions to any linear system of partial differential 
equations again to solutions provided the system is invariant under the group (as is the case for 
the time-harmonic Maxwell equations). Moreover it automatically yields solutions to the design 


equations with group G, because for any h E G, formula (147) implies via an elementary change of 


variables that hEo = X-a{h)Rio. Unfortunately, this simple method to obtain solutions to the design 
equations cannot be extended in a straightforward manner to non-compact groups which include 
translations. And projecting any special class of solutions to Maxwell’s equations is of course not 
guaranteed to deliver all solutions to the design equations. 
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15 Conclusions and outlook 


We have shown on the level of modelling and simnlation that twisted X-ray waves would be a very 
promising tool for structure analysis. 

Numerous theoretical challenges remain. For helical structures, a better understanding of the 
outgoing radiation in non-axial direction would be desirable. A related issue is to develop a general 
mathematical theory of the “radiation transform”, eq. (30), which takes the role of the Fourier trans¬ 
form when the incoming radiation is not given by plane waves. Robust phase retrieval algorithms 
need to be developed for structure reconstruction from twisted X-ray data. And for structures 
generated by non-abelian symmetry groups such as buckyballs, the right incoming waveforms are 
not clear; for reasons discussed at the end of Section in this case the design equations may not 
be the right approach. 

Challenges for the experimental realization include: generation of tunable coherent twisted X- 
ray waves with a broad spectrum of angular wavenumbers; axial alignment of incoming wave and 
structure; and achieving a sufficiently strong outgoing signal. 


16 Appendix: Intensity of outgoing radiation 

Here we relate the intensity ( [2^ of the real-valued diffracted electromagnetic radiation caused 
by any incoming time-harmonic held, given by the real part of eq. (0. to the complex electric 
held amplitude. We rely on the following simple lemma which is applicable due to the special, 
locally plane-wave-like features that the held vector 'Ejout is always perpendicular to k' and ^out is 
proportional to the vector product of k' and Eout- 


Lemma 16.1 (Poynting vector and intensity of locally plane-wave-like radiation). Sup¬ 
pose that 

E(x, f) = Re (E(x)e“*‘^*)), B(x, f) = Re (k(x) x E(x, f)) 

for some complex vector field E : —)■ and some real vector field k : —)■ with k(x) ■ 

E(x) = 0. Then 

S = —k f iRe EP cos^(a;f) -|- IlmEP sin^(a;f) -|- 2ReE ■ ImE cos(a;f) sin(a;f)^ (148) 

/io V / 


and 


I 




(149) 


Thus the Poynting vector of the real part of a locally plane-wave-like electromagnetic held is a 
huctuating mixture of contributions from both the real and imaginary parts of the spatial electric 
held, and the intensity is proportional to the absolute value squared of the complex electric held 
amplitude. One can hence say that the mystery of complex held amplitudes appearing in intensities 
of real, physical helds is explained by time-averaging. 


Proof Writing ReE = Ei, ImE = E 2 , we have E = Ei cos(a;f) — E 2 sin(a;f), B = (kx Ei) cos(a;f) — 
(k X E 2 )sin(a;t). Eq. (148) now follows from the expansion rule for double vector products, 
ax (b X c) = b(a ■ c) — c(a ■ b), and the orthogonality assumption k ■ Ei = k • E 2 = 0. Moreover 
since the time-dependent factor in (148) is always nonnegative, the absolute value |S| is a scalar 
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multiple of this factor and it is straightforward to evaluate the asymptotic time average in (22), 
which amounts to averaging over one period. Since cos^ and sin^ have average 1/2 and cos-sin has 
average 0, we obtain (149). 

The intensity of diffracted physical radiation at the observation point x, i.e. the intensity of the 
real part of the elecromagnetic field ( [TT] ), is therefore, by eq. (149), 

^ |Eo„i(x,f)|^ (eq. ([^, Section]^, 


where Eout is the complex electric held in (17) (note that its absolute value is independent of t). 
Here we have used the relations |k'| = cu/c and = \/(?. 
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